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Abstract 



Two consistency conditions for partition functions established by Akemann and Dam- 
gaard in their studies of the fermionic mass dependence of the QCD partition function 
at low energy (a la Leutwiller-Smilga-Verbaarschot) are interpreted in terms of integrable 
hierarchies. Their algebraic relation is shown to be a consequence of Wick's theorem for 
2d fermionic correlators (Hirota identities) in the special case of the 2-reductions of the 
KP hierarchy (that is KdV/mKdV). The consistency condition involving derivatives is 
• an incarnation of the string equation associated with the particular matrix model (the 

particular kind of the Kac-Schwarz operator). 



1. Recently G. Akemann and P.Damgaard Q established two consistency relations for the 
partition functions of matrix models which connect finite- volume partition functions with dif- 
i| ferent fermion numbers. In the course of their paper the similarity of these relations with 

different bilinear relations for the r-functions of the KP hierarchy were noted and they conjec- 
tured they were in fact related. This letter provides a comment on the origin of these reduction 
formulae. We show they are direct corollaries of Wick's theorem and the string equations of 
the r-function theory. Application of this theorem is possible because the partition functions 
of Cartanian matrix models are KP/Toda r-functions (see [|[ || and references therein, espe- 
cially § and §, for the generic theory and (7[ || for applications to unitary-matrix models). 
The algebraic relation ( "Consistency Condition I" of ref. jjj ) is true for any 2-reduction of KP 
r-function, which are essentially the KdV or mKdV reductions^, of which the Brezin-Gross— 
Witten (BGW) model || under actual examination in jjj is an example. The differential 
relation ( "Consistency Condition II" ) is the string equation, depending on the specifics of the 
particular model. (The BGW model is studied from the perspective of integrability theory in 
rcf.§.) 

These consistency conditions are imposed onto the QCD partition function as a function 
of fermionic masses /ij at very low energies, which is of the general form 

2& N '\{fi) = det<I>(M)/A({y})- (1) 
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1 The KdV and mKdV hierarchies are basically the same reduction of the (m)KP hierarchy. One deals with 
the mKdV case when one additionally considers the discrete (Toda) time "evolution". While generically the 
evolution of the KP r-function r n in this discrete time is governed by the Toda dynamics, in the case of the 
2-reduction there exist only two different values of the zero time due to the periodic boundary conditions, 
Tn+2 = r n . The corresponding r-functions tq and tj give the mKdV hierarchy. 



Here the matrix $ is defined by 

= <l>i(pi), h 3 = !>••• ,•#/) 



(2) 



for some functions <fo(/z) yet to be specified. The denominator is given by the Vandermonde 
determinant of the squared masses 
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A({ M 2 }) = Uitf-tf) = det^fr 1 ] 
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Now the consistency conditions I and II of jlj acquire the forms respectively 
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In what follows we discuss these equations, their origins and solutions. 

2. We begin with the universal (in the world of 2-reductions of KP) algebraic relation 
It states, in particular, that for certain sets of functions 0i(C)i i = 1, 2, . . . , and for any k 
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Here z, j = 1, . . . , 2k, while a, b = 1, . . . , k and the set of 2k variables {Q} is (arbitrarily) split 
into two subsets {£ } and {rib} of k variables each (for example, ( a = Ca+fc = 

Obviously, relation (^|) implies that all the functions </>i(C) are expressed in a certain way 
through the first two, </>i(C) and 02(C)- Sufficient conditions for (||) to be true are, for example, 
the recurrence relations: 



i+1 

j'=i 



const. 



(7) 



In the r-function theory this recursion corresponds to the (m)KdV reduction. The derivation 
of (^) from (Q) is an easy algebraic exercise (it can be done straightforwardly by induction 
in k or one can expand (f>i(Q — O -1 Sfc°=o-P ifc C _ ' c as a formal series in £ _1 and express the 
determinants through the characters of SL(N), see refs.|7| || for examples). 

Observe that, among other equations, the consistency condition I (||) and similarly (0) are 
invariant with respect to multiplication of all the 4>i(p) by an arbitrary function f(fi). 



3. Wc shall now comment on how (Q) and, more generally, (||), follow from the r- function 
theory. There are three ingredients needed from that theory: 

1. Wick's theorem and its realisation (the Hirota equation) in terms of the KP/Toda r- 
function in Miwa coordinates (involving positive and negative shifts of time-variables, 
see below); 

2. Determinant formulae for r in Miwa coordinates (with positive shifts only); 

3. The fact that r for a p- reduced hierarchy is essentially independent of the times t p k- 
Specifically, the (m)KdV r-function is independent of even times: 

dT KdV /dt 2n = 0. (8) 

Generally M, ^] the r-function is only defined up to times linear in the exponential: 
d 2 logr/dt^f, — 0. This is equivalent to the freedom to multiply the functions 4>i((J>) by a 
function f{p) mentioned earlier. Here we are making the specific choice that yields (|^). 

This last property allows one to essentially identify the positive and negative Miwa shifts of 
time-variables for the (m)KdV r-functions and combine Wick's theorem and the determinant 
representation into a common identity (0) as we shall now demonstrate. 



4. Wick's theorem (in the restricted form which we need here) is an identity between cor- 
relators for any free-fermion theory with action of the form J 7 ?)V'( 7 7)^' D £ C ^ D? 7- It 
says 



n>te>)fc)> = det^^a)^)). (9) 

o=l 



For the KP/Toda r-functions, when D = 2 and G(£,r)) is a first-order differential operator, 
C(Cj r ?) = ^(£7 + ^4(0): the correlators can be written as 

*Ji T K p{t) n a ,6(£«-w 

The dependence on times {t} on the left-hand-side of this expression is hidden in the definition 
of the brackets. Here the arguments of the r-function are Miwa-shifted time- variables 



and the Vandermonde determinants (A(£) = Ila<b(£a — &>) e ^c.) arise from the normal 
orderings. While 

tkp (t + - foj) j = r KP (t)( : Y[ ^(LMVa) :) 

is non-singular at coincident points £ and 77, the correlator in ( |l0| ) vanishes when any two of 
£ or any two of rj collide and have a pole when any £ approaches any rj. The Vandermonde 
determinants encode this singular behaviour and corresponds to the naive ( "Riemann sphere" ) 
normal ordering. One sometimes uses more sophisticated orderings associated with non-trivial 
Riemann surfaces, but this is not the case (for the matrix models) under consideration. 



5. In Wick's theorem, the numbers of £ and r\ variables (i.e. of the ip an d tp operators, or 
of the positive and negative Miwa shifts) are necessarily the same. However, one can take a 
limit when all the ij points collide, say at infinity. Then Wick's theorem implies a peculiar 
determinant formula for the r-function in Miwa parametrisation || (with only positive Miwa 
shifts remaining): 

T KP (t) A(0 ' 1 ' 

From this perspective, the functions (j>i(£) arise as expansion coefficients in 

oo 
i=0 

Note that all the time-dependence is hidden in the shapes of the functions fy. In particular, 
one can make a further Miwa transformation of the times t in ( |l^ ) enlarging the right-hand 
determinant of ( |i1^ ) and correspondingly changing the functions The same is true of the 
//-dependency in formulae (^)-@: one could hide the dependence on {/t.;} in the functions (f>j, 
though they can of course be taken out to produce determinant formulae with additional sets 
of Nf parameters fit,. ■ . , Hn } ■ 

At this stage nothing is gained by substituting the determinant formula ( |l2|) back into the 
Wick identity (^|), from which it was derived: expansion of both sides of (j^) in powers of ?/~ 4a 
will produce nothing new. 



6. However, in the case of (m)KdV r-functions there is an amusing way of introducing (an- 
other set of) r\ parameters back into the game. Namely, because of property (||), 

T Kd v{t + [£] - [-»?]) = r KdV (t + [£] + [t]}) (13) 



since the difference between the arguments is seen (by (|llj)) only in the even times t2n, on 
which TKdv does not depend. Now one can apply Wick's theorem (S) with the substitution 
( [l0| ) to the left-hand-side of ( |l3| ) and the determinant formula (|12) to its right-hand-side. 
Straightforward calculation produces the relation (^j) . 

Note now that for the ratio of determinants ( |l2| ) to be a r-function, the functions 4n (ji) are 
to have prescribed asymptotics at large /i, usually <fii{p) ~ (see, however, [Q). This en- 
sures that the r- function is an (iVy-independent) function of the usual times f& = 1/k Y]- s , 
and so allows an expansion as a formal series in the time variables tj~ . These asymptotic prop- 
erties remain unchanged when the functions 4>i(p) are multiplied by a function /(/it) that 
behaves like 1 for large fi. In particular, if log/(/x) expands into a Taylor series in such a 
redefined r-function is multiplied by a linear exponential of times t/-. Therefore, the property 
(0) holds for the (m)KdV r-function only for a very specific normalisation of the 4>i{^). 

We have now shown how consistency condition I (0) can be understood in terms of Wick's 
theorem for a particular (m)KdV r-function. It remains to discuss consistency condition II 
and the matrix theory origins. 



7. While the algebraic relation (|J) allows one to impose restrictions on <pi (//) that lead to an 
(m)KdV r-function, i.e. express all ^(/i)'s through the first two, the differential relation (||) 
expresses all of the </>i(/z)'s in terms of derivatives of the first. The combination of these two 
restrictions then fixes the 4>i(p) (up to a linear combination of rows in the determinant and a 
trivial prefactor). In terms of integrable hierarchies, this relation is called the string equation 



|1 . The role of this equation is to fix a unique solution of the integrable hierarchy (specified 



by the algebraic relation (|4|) and the form of Z, 

Let us start with the simplest case of Nf = 0. Then, '(§) implies that|0 2 (/x) ~ fidcf)i(fj,)/dfj,. 
This is, however, still not a restriction on 4>i(n) itself. In order to see how (||) restricts all 
the <^i(/i), one needs to consider it at all values of Nf. In doing so, it is necessary to use the 
equation that is tautologically satisfied by the ratio of determinants ( p^ ) [||[]: 

T Nf+2 ({fi},£,v) = (£-r])T N .({[j,}) [ TAr /+ 1 ^ /i J''^ fjv /+ 1 ^ /i J''^ ~~ ^ ^ ^] ^ 14 ) 

where the hat over tn means that the last row <^Af({/i}) is substituted with </>jv+i({/i}). This 
equation is typically used for derivation of the string equation EL |^] . 

Now, since J2i ^i^^i^N f (p 2 ) — Nf(Nf — l)Ajv.(// 2 ), one can rewrite (Bh as 



det 
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Comparing with (|l4|) , we conclude that 



+ ^ ) det^-t- 1 ^^}, 



det^+^M,??)-^^) 



(15) 



det (JV ' +1) l>(M,£)det (A, ' +1) $(M,77) -(£<-► ??). 



with the evident notation $ denoting the matrix whose last row </>jv({/i}) has been substituted 
with cj)N + i({fi}). Induction on Nf (we have already shown it true for Nf — above) shows 
this equation is solved if 



detl> = ( E j det $ + N f a det$, 



(16) 



a being an arbitrary constant. In fact this relation is another version of the string equation 
(cf. with ||). The term proportional to det$ can always be absorbed by the simultaneous 
rescaling all <fii{p) with /i Q and ( |l6|) yields in its turn that 4>i(n) ~ (/xc^) 1-1 ^^). 

Therefore, the consistency condition II really restricts the functions 4>i(n). Combining the 
form 4>i(ii) ~ (/i9 A1 ) 4_1 0i(/Lt) together with the reduction condition (]?]) (^consistency condition 
I), one arrives at the solution M 



4>i(n) ~ fi a+l - 1 I 1/+i _ 1 (fj,) 



(17) 



where Ik(p) are the modified Bessel functions, v and a are arbitrary parameters (our solution 
slightly differs from that in JjJ by the additional parameter a) . 

2 Up to a term proportional to <j>i{i±)\ there is always a freedom to add to any ipk( x ) a line ar combination 
of <t>i(^i) with i < k with constant coefficients, since it does not change the determinant in (|f 2[). Hereafter, we 
denote this freedom with the sign 

3 The authors of 111 correctly suggested its relevance to the problem. However, let us emphasise again that 
this equation does not select out any particular r-function, or any particular 4>i(fJ.). This is just another form 
of formula (Q. 



8. The solution ( |l7| ) obtained for can also be presented in terms of the matrix integral 

I 

zl Nt \{n})~ [ dXe T ^ M2x/ ^ {l/+Nf)loeX+1/x+{a ^ l ' )losM '\ (18) 



where the integral is taken over Hermitian Nf x Nf matrix X and \ii are the eigenvalues of the 
matrix M. For the particular values a = v = this integral may be rewritten as an integral 
over unitary Nf x Nf matrix U [0 

Z$, Nf) m)~ [dUe Tr [ J ' u+uij \ M 2 = J,P. (19) 



The function Z„ f ({jti}) by itself is certainly not a r-function for the Miwa coordinates 
(11) we are using here (chosen because of the formulation of [1]). As we noted earlier, for the 
ratio of determinants (O) to be a (r-)function of the times tk we required the asymptotics of 
4>i(n) at large /i to be normalised by ||, which are not those of (|l7|). An easy calculation 
shows that, in order to obtain a r-function, one needs to include several prefactors j^j 

^(W)=II(^ + ^)n^Mr a e-' li ^ ) ({M}). (20) 

i<j i 

Here the parameter v plays role of the discrete time (in the standard notation of j| it has 
the opposite sign). Such prefactors typically emerge in matrix model theory: to obtain a 
r-function, one should factor out a quasi-classical piece from the matrix integral. For the 
particular integral ( |l8| ) this exactly gives (|2^), and in fact reference Q also shows the inde- 
pendence of even times we have wanted here. Thus we return to the matrix model integrals 
that inspired the consistency conditions we have been discussing. 

In concluding let us observe that these integrals can be further generalised yielding other 
r-functions. In particular, one can consider an arbitrary function V{X) polynomial in X or 
1/X. Then, the integral 

C V J dXe ^ XL -^ +N ^ x+v ^ (21) 
results in a two-dimensional Toda lattice r-function, 



dcty 



r v {{t}) = L_ 1 (22) 

Here CV is the normalisation factor cancelling the quasi-classical part of the integral, and \ii 
are related with the eigenvalues k of the matrix L via the saddle point condition I + V'(p) = 
or l + V'(l/fj,) = depending on whether V(X) is a polynomial of X or 1/X respectively. The 
time variables are again tk = /i^ fe and 

= 1 e vW-*v'W f dxx W) e *i+vW (23 ) 
y/V"((J.) J 

where (i, (3) is equal to i — j3 — 1 for V(X) being a polynomial in X and (3 — i for a polynomial 
in 1/X. 
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